Electron-phonon coupling in the conventional superconductor YNi 2 B 2 C at high 
phonon energies studied by time-of-flight neutron spectroscopy 
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We report an inelastic neutron scattering investigation of phonons with energies up to 159 meV 
in the conventional superconductor YNi2B2C. Using the SWEEP mode, a newly developed time- 
of-flight technique involving the continuous rotation of a single crystal specimen, allowed us to 
measure a four dimensional volume in (Q, E) space and, thus, determine the dispersion surface and 
linewidths of the A\ g (w 102 meV) and A u (« 159 meV) type phonon modes for the whole Brillouin 
zone. Despite of having linewidths of V = lOmeV, Ai g modes do not strongly contribute to the 
total electron-phonon coupling constant A. However, experimental linewidths show a remarkable 
agreement with ab-initio calculations over the complete phonon energy range demonstrating the 
accuracy of such calculations in a rare comparison to a comprehensive experimental data set. 

PACS numbers: 74.25.Kc, 78.70.Nx, 63.20.kd, 63.20.dk 



Superconductivity is one of the most intriguing phe- 
nomena in condensed matter physics and each discovery 
of new compounds, such as the iron based superconduc- 
tors, sparks a tremendous amount of scientific interest. 
To date, BCS theory Q, with the Cooper pairing mech- 
anism provided by electron-phonon interaction, repre- 
sents the only well-understood and experimentally ver- 
ified microscopic picture of superconductivity, and sys- 
tems which are successfully described by BCS theory 
are known as conventional superconductors. The cor- 
responding superconducting transition temperature can 
be calculated by, e.g., the Allen-Dynes formula 0]: 



1.2 



exp 



-1.04(1 + A) 
A-/x*(l + 0.62A) 



(1) 



Here, A is the electron-phonon coupling (EPC) constant, 
H* the effective electron-electron interaction and u e f f 
the effective phonon frequency. Thus, a large effective 
phonon frequency and electron-phonon coupling constant 
are favorable for a high superconducting transition tem- 
perature. Determining A and w e ff microscopically, i.e. 
by measuring all phonon energies and linewidths, is ex- 
tremely time consuming if not impossible. Therefore, 
ab-initio calculations of the lattice dynamical properties 
are often used to estimate A and the corresponding su- 
perconducting transition temperature T c [3]. 

Following the discovery of superconductivity in the 
rare-earth nickel-borocarbides (RNBC) the conven- 
tional phonon mediated pairing mechanism was quickly 
accepted to be valid for these compounds, based on 
density-functional-theory (DFT) predicting a large A = 
2.6 [H and the observation of a strong isotope effect 



The calculated electronic structure [B| showed that Ni 
and B states dominate at the Fermi energy. When a 
strong sensitivity of T c towards the bonding angle of the 
MB4 tetrahedron was observed [8| , it was even suggested 
that high energy phonons of A\ g symmetry, which mod- 
ulate this bonding angle, were prominently involved in 
mediating superconductivity This is somewhat 

surprising, as, according to equation [TJ an EPC constant 
of A = 1 would only require uj e f f = 20 meV in order to 
explain a T c of about 15 K (fi* = 0.13), as it is observed 
in the RNBC compound YM2B2C. In contrast, the fre- 
quency of the Ai g mode is close to 100 meV [lfjj. 

In this letter, we report measurements of the Ai g and 
A u type phonons at very high energies, i.e. E w 100 meV 
and 160 meV, respectively, and another pure B-C mode 
(see supplemental information) dispersing at 37meV < 
E < 52meV in YM2B2C. Utilizing a novel technique 
whereby the sample is continuously rotated, the SWEEP 
mode developed at the TOF chopper spectrometer ARCS 
[ill ] located at the Spallation Neutron Source (SNS), Oak 
Ridge National Laboratory, we were able to sample a four 
dimensional volume in (Q, E) space, from which we ex- 
tracted constant energy planes in reciprocal space (Fig.[T]) 
or data for varying energy transfers and constant Q along 
all high- and off-symmetry directions (Fig. [2]). Results 
from TOF spectroscopy are complemented by triple-axis 
measurements for lower energy phonon linewidths requir- 
ing an improved energy resolution. We compare our 
results to detailed calculations of the lattice dynamical 
properties based on density-functional-perturbation the- 
ory (DFPT), which recently challenged the importance 
of the high energy modes for superconductivity in boron 



2 



o 

x~ 
x~ 




10 15 
(0, 0, L) (r.l.u.) 



20 



FIG. 1. (Color online) Raw data collected in the (hhO) - (001) 
plane with an incident energy of Ei = 160 me V. The figure 
shows a cut at E = (102.5±2.5) meV. The data were obtained 
with an angular rotation of the specimen of 60°. The black 
pixels illustrate the wave vector range covered with an angular 
range of 1° (no intensity scale applied). Note that the wave 
vectors covered in an experiment with such a small angular 
range changes significantly with energy transfer. 



carbides [12J. 

Neutron scattering experiments for such high energy 
excitations are challenging as one has to measure a sig- 
nificant part of the Brillouin zone with an energy reso- 
lution comparable to the electronic contribution to the 
phonon linewidth. The low neutron flux at energy trans- 
fers of about 100 meV and concomitant poor energy res- 
olution exclude the usage of thermal neutron scattering 
on a triple-axis spectrometer (TAS), which had nearly a 
monopoly on energy and momentum resolved phonon in- 
vestigations up to very recent times. In the last decade, 
high-resolution inelastic x-ray scattering (IXS) has be- 
come a serious competitor to TAS measurements. How- 
ever, phonons at high energy transfers are generally as- 
sociated with displacements of the light elements, which 
have only small x-ray scattering amplitudes. 

Time-of-flight (TOF) chopper spectrometers at mod- 
ern pulsed neutron sources [13| offer a way out of this 
predicament. They combine a high neutron flux over a 
large energy range with good energy resolution at large 
energy transfers 



11] 



Furthermore, these instruments 
map large areas of reciprocal space and energy simul- 
taneously employing arrays of position sensitive 3 He de- 
tector tubes of up to 30 m 2 ll|, |l4|. However, the layout 
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FIG. 2. (Color online) (a)(c) Color-coded contour plot of raw 
data for varying energy transfers vs wave vectors Q = (0, 0, I) 
obtained with incident energies of Ei — 160 meV (left) and 
250 meV (right). (b)(d) Examples of data analysis of energy 
scans extracted from the measured (Q, E) space at the given 
wave vectors and incident energies of Ei = 160 meV (left) and 
250 meV (right). Constant-Q scans are binned over ±0.1 r.l.u. 
along (00L) and ±0.2 r.l.u. along (HHO) and (HHO) and are 
analyzed by fitting a Lorentzian convoluted with the instru- 
mental resolution function. Solid/dashed/dotted lines are the 
respective convoluted fit/resolution function/experimental 
background. Experimentally observed phonon frequencies 
and linewidths (FWHM) of the Ai 9 (filled symbols) and A u 
(open symbols) phonon modes are shown in panels (e) and 
(f). Solid and dashed lines were calculated from DFPT. The 
dotted part of the line in (f) indicates that these values could 
be calculated only approximately (see text). 



covering a cylindrical surface around the sample posi- 
tion does not allow one to make energy scans at constant 
momentum transfer Q, the usual choice in neutron TAS 
measurements. Instead, the instrument samples a curved 
surface of reciprocal space as indicated in Fig. [TJ which 
makes it impossible to directly extract a lifetime of an ob- 
served excitation. To do so requires fitting of a full four- 
dimensional resolution convolution with an a-priori as- 
sumed model for the dispersion relation . In addition, 
the range of measured wave vectors strongly depends 
on the energy transfer. In measurements of quasi two- 
dimensional, e.g., magnetic excitations, this drawback is 
not too serious because all results may be projected onto 
a basal plane. For phonon measurements however, and 
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generally for measuring any excitation that is strongly 
dispersing in all three dimensions, such a projection is 
rarely if ever meaningful. This has stimulated the de- 
velopment of a new technique for the determination of 
phonon dispersion curves using a TOF instrument, which 
we call the SWEEP mode. The SWEEP mode involves 
a continuous rotation of the single crystal sample dur- 
ing the measurements so that large volumes of reciprocal 
space are sampled. We note that the SWEEP mode pro- 
vides greater flexibility than the related HORACE mode 
currently being developed at ISIS [l6| , where the sample 
is measured at a series of discrete angular positions over 
a period of a day or more. In the SWEEP mode, the 
entire volume is measured within minutes, allowing the 
experimenter to optimize the angular range as the data 
are accumulated. 

The sample we used was the same single crystal as 
used previously in the TAS measurements of the low en- 
ergy phonons P3, EH with a volume of about 0.4 cm 3 . 
Thus, all the progress made in the TOF measurements 
came from the technique, and not an increase in sam- 
ple mass. Data were collected using incident energies of 
Ei = 160 meV and 250 meV for measuring the branches 
involving the A lg (ss 100 meV) and the A u (« 160 meV) 
modes, respectively. Because these modes are flat in mo- 
mentum space, their lifetimes can be easily extracted 
from energy cuts at constant Q, obtained from the re- 
construction of the full S(Q, E) 1 and only require a one- 
dimensional resolution convolution. The resolution func- 
tion, an Ikeda-Carpentcr function describing the moder- 
ator [i(3] convoluted with a Gaussian describing the chop- 
per contribution, was determined at zero energy transfer 
from vanadium scans and scaled to finite energy transfers 
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The so determined instrumental full widths at half 
maximum are 4.0 meV (Ei = 160 meV, AE = 102 meV) 
and 7.4 meV (E t = 250meV,AE = 159 meV) (Figs. Eb, 
d). The sample was mounted in a closed-cycle refrig- 
erator and measurements were done at T = 5 K. The 
total data acquisition time for the results presented here 
were 6 days. The sample was continuously rotated back 
and forth throughout the experiment with an angu- 
lar speed of about 5° per minute over a range of 60° 
(Ei = 160 meV) or 38° {E t = 250 meV). Each de- 
tected neutron was stored with the respective sample 
rotation angle in order to allow a later transformation 
into (Q, E) space. Linewidths of the mode dispersing 
at 37meV < E < 52meV were obtained from standard 
triple-axis experiments on the IT and PUMA spectrome- 
ters located at the ORPHEE reactor at LLB, Saclay, and 
at the research reactor FRM II in Munich, respectively 
(for experimental details see supplemental information). 

Raw data obtained in TOF spectroscopy are shown 
in Figs. [T] and [2] Obviously, there is a strong variation 
of the scattered intensities for both of the high energy 
phonon branches as a function of wave vector along the 
line (00.L), whereas there is nearly no variation along 
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FIG. 3. (Color online) Left-hand scale: Ab-initio calculations 
of the phonon density-of-states (black, solid line), the Eliash- 
berg function a 2 F(uj) (blue, dashed line). Right-hand scale: 
Partially energy integrated electron-phonon coupling constant 

A = 2 ■ duj 1 (red, solid line). 



(HH0) (Figs. □ and UK, c). We note that the observed 
intensity variations quantitatively agree with that pre- 
dicted by DFPT based on the calculated phonon pat- 
terns (see supplemental information). For a detailed 
comparison of the observed and calculated frequencies 
and linewidths, constant Q scans, binned over ±0.1r.l.u. 
along (00L) and ±0.2r.l.u. along (HH0) and (HH0), 
were fitted by a Lorentzian convoluted with the instru- 
mental resolution function [wave vectors are given in re- 
ciprocal lattice units (r.l.u.) of (2n/a, 2n/b, 2n/c), where 
a = b = 3.51 A, and c = 10.53 A]. Two examples of 
fits to the raw data are shown in Figs. [3]Jb) and (d). 
In this way, we investigated large sections of the Bril- 
louin zone including the main high symmetry directions 
along the crystallographic (100), (110) and (001) axes. 
Representatively, we show results for the (001) direction 
in Figs. EJe) and (f). For both phonon branches, we 
see good agreement between calculated and observed fre- 
quencies with deviations of less than 2% over the whole 
Brillouin zone. What is more, the calculated and ob- 
served phonon linewidths, which are all in the range of 
8—10 meV, agree quite well for the branch starting from 
the Ai g mode. Hence, there is strong EPC throughout 
the Brillouin zone. Why it nonetheless does not signif- 
icantly contribute to A will be discussed below. As to 
the branch starting from the A u -mode, theory and ex- 
periment agree well except for the (001) direction: the 
phonon is resolution limited in most parts of the Brillouin 
zone, but acquires a detectable linewidth near the zone 
boundary along the (001) direction. This is qualitatively 
predicted by theory, but quantitatively DFPT overesti- 
mates the phonon linewidths, i.e. DFPT will predict a 
larger impact of the A u modes EPC than experimentally 
corroborated. We now use the DFPT results for a discus- 
sion about the contribution of the various phonon modes 
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FIG. 4. (Color online) Observed (a) and calculated (b) disper- 
sion surface of a phonon mode showing a pronounced anomaly 
in the (hhO)- (001) plane. Experimental phonon energies were 
determined from TOF data in steps of 0.1 r.l.u. in h and 
/. (c)(d) Measured and calculated phonon linewidths of the 
mode shown in (a) and (b), respectively, along two high sym- 
metry lines of the Brillouin zone. Experimental results shown 
in (c) and (d) were obtained by triple-axis spectroscopy (see 
text). 



to the total electron-phonon coupling constant A. 

In spite of the large phonon linewidth of the A\ g and, 
to some extend, A u modes, DFPT claims that two thirds 
of the total A is due to phonons with E < 30meV 
(Fig. 13]) and in fact more than 97% is due to phonons 
with E < 70 meV. In order to understand why the size- 
able linewidth of the A\ g mode does not contribute more 
strongly to A and thereby to superconductivity, it is in- 
structive to look into the different components of the 
Allen-Dynes formula (equation [T]) . From the definitions 
of the EPC coupling constant A, 



A = 2 • / duj 



and the effective phonon frequency w e //, 



(2) 



J eff 



exp\-- 



a 2 F(cu) . 



ln(u) , (3) 



we see that both depend on the Eliashberg function 
a 2 F(uj) divided by the frequency. The importance of 
the weighting factor 1/uj can be seen in Fig. 02 which 
shows the calculated phonon density-of-states (PDOS), 
a 2 F(co) and also the partially energy integrated electron- 
phonon coupling constant A. The A\ g modes around 



100 meV clearly have an enlarged contribution to a 2 F(ui) 
with respect to their spectral weight in the PDOS, in- 
dicative of substantial EPC. But in the end, their in- 
tegrated contribution to A remains fairly small, much 
smaller than that of phonon modes below 25 meV: there- 
fore, phonons above 70 meV account only for 3% of the 
total A. Likewise, the value of w e ff is only slightly en- 
hanced by the high-frequency modes: the value computed 
from the DFPT results is w e// = 23.16 meV. There- 
fore, these modes are not decisive for the relatively high 



15.6 K. These conclusions are supported by 



the above mentioned fact that the observed linewidths 
of the ^4 u -related modes are even smaller than predicted 
by DFPT. 

While DFPT clearly shows the negligible role of the 
high energy phonon branches and correct ly pre dicts 
strong EPC in the acoustic phonon branches [17l.ll8ll22|. 
it also predicts an extremely steep phonon dispersion of 
a pure Ni-B mode dispersing from 52 meV at the zone 
center to 36meV at the Z point, i.e. q = (0,0, 1). This 
anomalous dispersion, with a minimum centered at the 
Z point, is clearly related to strong EPC (see Fig. @] and 
Ref. [T2I ) . Experimentally, we were able to extract corre- 
sponding frequencies over a large section of the Brillouin 
zone from our data set with Ei = 160 meV. We compare 
the obtained dispersion surface in the (hh0)-(00l) scatter- 
ing plane, consisting of 169 individually determined fre- 
quencies, to the calculated dispersion in Fig. fU We note 
that it would be very time consuming and impractical 
to use a conventional TAS for the same purpose even at 
these lower energy transfers. As can be seen from a com- 
parison of Figs.QJa) and (b), the experiment and theory 
agree quite well. Due to the limited amount of beam time 
we were not able to investigate the corresponding phonon 
linewidths requiring better energy resolution, i.e. a lower 
incident energy. Therefore, we show phonon linewidths 
along (h, 0,0) and (h,h,l) for the mode dispersing be- 
tween 37 and 52 meV obtained in a TAS measurement in 
Fig. ESfcJ and (d)[1 ]. We observe good agreement be- 
tween the predicted contribution to the phonon linewidth 
and the experimental results. Taking into account pre- 
dicted and reported strong EPC in acoustic branches of 
YNI2B2C [la ], our new results for phonon modes with 
37meV < E < 160 meV demonstrate that DFPT yields 
reliable predictions of the lattice dynamical properties 
over the whole energy and wave vector spectrum in the 
conventional superconductor YNi 2 B 2 C. 

In conclusion, the recently developed SWEEP mode in- 
volving a continuous rotation of a single crystalline sam- 
ple on a TOF chopper spectrometer was used to obtain 
very detailed information on phonon frequencies as well 
as phonon linewidths over a large region in reciprocal 
space and up to very high phonon energies using a sam- 
ple of moderate volume. Using this technique, we show 
that the high energy vibrations in the superconductor 
YM2B2C couple indeed very strongly to the electrons, 
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but that this coupling is nevertheless irrelevant to su- 
perconductivity, in contrast to earlier claims 0, Q and 
the common association that high energy phonons with 
strong electron-phonon coupling yield a large effective 
phonon frequency ui e ff and, consequently, T c . The latter 
only applies for superconductors like MgB2, where nearly 
all of the EPC is in one high energy branch [23], but it 
does not hold true if there are several phonons, both at 
low and high energies, with strong EPC, as is the case 
for YN12B2C. Here, as presented in detail above, the ef- 
fective phonon frequency is dominated by the low energy 
phonons, and the high energy phonons provide only a 
small contribution. Our results show good agreement 
with DFPT calculations for the phonon energies and 
linewidths and corroborate the DFPT result that lower 
energy phonon modes with large amplitudes of the lighter 
atoms mainly mediate superconductivity in YN12B2C. 
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